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Introduction: The issue of whether crustal con- 
tamination has affected the lithophile trace element 
budget of shergottites has been a point of contention 
for decades [1-7]. The evaluation has focused on the 
enriched shergottite compositions as an outcome of 
crustal contamination of mantle-derived parent mag- 
mas or, alternatively, the compositions of these stones 
reflect an incompatible trace element (ITE) enriched 
mantle source. 

Recent findings have supported the notion that en- 
riched shergottite compositions reflect crustal contam- 
ination. For example, ITE-enriched olivine-hosted melt 
inclusions have been found in depleted shergottites [7], 
the rare earth element (REE) compositions of Mars 
crustal sample NWA 7034 and its paired stones, when 
mixed with depleted shergottite compositions, ade- 
quately reproduce the bulk REE compositions of en- 
riched shergottites [8,9]. These two results are con- 
sistent with findings from Cl and S isotopes that indi- 
cate shergottite magmas interact with martian crust on 
their way to or, on the surface [10,11]. However, a 
crustal contamination model must adequately explain 
not only the ITE compositions of enriched shergottites, 
but also their lithophile element isotope compositions 
such as for Nd and Sr isotopes. Using Os-Nd isotope 
relationships, realistic crustal contamination models 
have been shown to not work to produce enriched 
shergottites from depleted shergottite parent magmas 
[12,13]. Hence the origin of the ITE-enrichment in 
enriched shergottites remains an open question. 

To further evaluate this issue, a comprehensive 
study on enriched shergottite NWA 856 [14] was un- 
dertaken, building on a previous ®’Rb-*’Sr and '47Sm- 
'3Nd isotope chronology study [15]. The new work 
included major and trace element analyses of minerals, 
and evaluated the crystallization history of NWA 856. 
This stone has a bulk molar Mg# of 0.49 and lacks 
olivine, making it one of the most evolved members of 
the enriched shergottites. Hence, this stone underwent 
extensive crystallization from its parent magma such 
that there was a strong prospect for crustal contamina- 
tion during its evolution. In addition, Nd isotope com- 
positions have recently been reported for Mars crustal 
breccia NWA 7034 [8,16]. Combining these new data 
with the comprehensive data available for shergottites 
allow for a refined evaluation of the effects of crustal 
contamination and its implications to shergottite source 
variation in the martian mantle and crust. 


Results: The previously reported *’Rb-*’Sr and 
47§m-!'3Nd results were re-evaluated for ages using 
ISOPLOT. The new computations give *’Rb-°’Sr and 
'47§m-'3Nd internal isochron ages of 162414 Ma (+20 
and used herein) and 162.7+5.5 Ma, respectively, with 
an initial eNd = -6.6+0.2. These ages fall in the range 
for those of the other enriched shergottites that when 
combined give an average Nd isotopic age of 173422 
Ma, equivalent to an average age of 175 Ma proposed 
by [6]. 

The bulk rock composition of NWA 856 was used 
in MELTS simulations to model equilibrium and frac- 
tional crystallization sequences to compare with the 
crystallization sequence from textural observations and 
to the mineral compositions. These models constrain 
the depth of initial crystallization to a pressure range of 
0.4-0.5 GPa (equivalent to 34-42 km) in anhydrous 
conditions at the Fayalite-Magnetite-Quartz buffer, 
and consistently reproduce the observed mineralogy 
throughout the sequence with progressive crystalliza- 
tion. The REE abundances in the clinopyroxenes and 
maskelynite are consistent with progressive crystalliza- 
tion in a closed system. Closed system crystallization 
also explains best the compositions of other enriched 
shergottites including LAR 06319 [17] and NWA 
5298 [18]. 

The Ti/Al ratios in the NWA 856 clinopyroxenes 
are consistent with initial crystallization occurring at 
the lower crust to uppermost mantle depths constrained 
by the MELTS models followed by polybaric crystalli- 
zation as the parent magma ascended to the surface. 
Polybaric crystallization is recorded in pyroxene Ti/Al 
ratios for a wide range of depleted and enriched sher- 
gottites including for example, depleted shergottite 
NWA 5789 [19], and enriched shergottites RBT 
04262, NWA 5398, LAR 06319, Zagami, Los Ange- 
les, and Shergotty (not shown). 

Discussion: The mechanisms of crustal contami- 
nation in terrestrial mafic magmas are well understood 
[20-22]. Crystallization provides latent heat that can 
then be used to melt the surrounding crust. Two end- 
member scenarios of crustal contamination are possi- 
ble, one where the compositional outcomes will follow 
trends of assimilation and contemporaneous fractional 
crystallization (AFC), or they will simply follow mix- 
ing trends between compositional end-members. Any 
successful model that ascribes compositional variabil- 


ity to any given related 
suite of mafic srocks to 
crustal contamination 
should thus be able to 
predict the outcomes for 
major element, trace 
element, and _ isotopic 
variability. As noted 
eee) ber above, the bulk and 
mineral compositions of 
depleted and enriched 
shergottites, including 
NWA 856, indicate that 
polybaric crystallization 
occurred while — sher- 
gottite parent magmas 
were transported 
through the martian 
crust and underwent variable amounts of differentia- 
tion. Hence there was ample opportunity for sher- 
gottite magmas to interact with the martian crust and 
inherit its compositional signatures either through AFC 
processes or mixing between shergottite parent mag- 
mas and crust. 

The trends for mixing between depleted shergotite 
compositions with martian crust represented by NWA 
7034 are compared to the compositional variation dis- 
played by all shergottites (Figure 1). For Ce versus 
Ce/Yb, both mixing and AFC models can explain the 
bulk rock variation observed in intermediate and en- 
riched shergottites as interaction of depleted sher- 
gottite magmas with martian crust. This is consistent 
with previous proposals that their REE budgets result- 
ed from crustal contamination [8,9]. The crust- 
depleted shergottite mixing and AFC models plot 
through the high molar Mg#-eNd data for enriched 
shergottites such as those for RBT 04262. However, 
the models fail to produce, for example, the restricted 
range of eNd at 150 Ma of -6.8+0.6 over a wide range 
of molar Mg# from 0.65 to 0.25 (Figure 1), or the 
range of Ce/Yb ratios versus molar Mg# (not shown). 
As the amount of crust proportions increases during 
progressive mixing or AFC in a differentiating sher- 
gottite magma, all of the models predict that the re- 
spective molar Mg# decreases, and that eNd would 
instead show a progressive decrease in the evolved 
magmas. Hence, the constant eNd of enriched sher- 
gottites over a wide range of differentiation indices is 
not an outcome of crustal contamination of a depleted 
shergottite parent magma. 

The only crustal contamination model that can ex- 
plain the bulk lithophile ITE compositions of enriched 
shergottites is one where primitive, high molar Mg# 
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Figure 1. AFC and mixing 
of crust with depleted sher- 
gottite models. Bulk rock 


magmas experience crustal contamination, that is then 
subsequently shut down or strongly minimized during 
the remaining differentiation. While not impossible, 
evidence from terrestrial systems show instead that 
crustal contamination is a progressive process during 
differentiation of mafic magmas [20-22]. Hence, it is 
concluded here that the bulk lithophile ITE budget of 
enriched shergottites represents instead a mantle 
source characteristic, likely formed during the early 
differentiation history of Mars [4-6]. 
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